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Abstract 

In the color-flavor-locked (CFL) phase of the QCD superconductor we show that 



> 

in 

m _ 

VO \ baryons behave as qualitons (called "superqualitons" ) with quantum numbers B = 

(1 mod 2)/3, S — 1/2 and Y = B. An intriguing possibility implied by this identi- 

Q\ ' fication is that light baryonic modes in the form of superqualitons could be excited 

below the (color) superconducting gap in the CFL phase, a novel phenomenon foreign 

Mh, to normal superconductors. 
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1. Dense hadronic matter poses still a challenging problem to theoretical physicists 
despite years of vested efforts. With the advent of dedicated fixed target accelerators 
(SIS) and heavy-ion colliders (SPS, RHIC, LHC) the theoretical issues have become more 
imperative as experiments may shed new light and spur new interest in the problem. 

At subnuclear matter densities, hadronic matter behaves as a free gas of nucleons 
and a good understanding of this phase can be reached by using virial type analyses |Q. A 
liquid-gas or Ising transition has been proposed |2J and recently observed ||] . At a few times 
nuclear matter densities, phenomenological analyses guided by nuclear structure data and 
constrained by neutron star observations have revealed a variety of phenomena including the 
possibility of a kaon-condensed phase that would be at the origin of proton-rich stars pj. 
At densities larger than nuclear matter densities it has been recently suggested that a 
variety of diquark-rich and superconducting phases may take place with potentially large 
gaps ||, y, 0, ||, || . Present neutron star cooling scenarios however appear to disfavor large 
gaps at moderate densities, although this conclusion is by no means definitive [1C]. 



At asymptotic nuclear matter densities, quarks interact weakly [|ll]] due to asymptotic 



freedom [12]. The ensuing metallic phase is supposedly screened except for possible infrared 



problems. It was originally suggested by Bailin and Love [ 13 1 and subsequently restressed by 
others [|| that in the asymptotic regime quarks pair at the Fermi surface resulting in small 
energy gaps. The perturbative arguments were recently revisited in light of the fact that the 
magnetic modes undergo Landau damping instead of screening, with the possibility of large 
gaps [||, [7|, H|, HJ. Here we just note that perturbation theory may still fail at asymptotic 
densities due to the nearness of the Gribov horizon in gauge theories [14]. 

Assuming that the superconducting phase is favored at asymptotic densities, it was 
recently observed that for three massless flavors, quarks may pair into diquarks with color 



and flavor locked forming the CFL (color-flavor-locked) phase [15]. It was further suggested 



that the CFL phase exhibits a spectrum with matching quantum numbers to that of the 



zero density phase [ 16 ] , a situation reminiscent of the strongly coupled standard model |lj 
and more intriguingly of the Cheshire-Cat phenomenon in low-energy hadronic structure 
(see [jl^] and references therein). The CFL phase appears to be energetically favored by 
effective theory calculations ||. 

In this letter we show that in the superconducting CFL phase, baryons behave as 
qualitons (called in short "superqualitons"), realizing Kaplan's scenario in the high density 
phase pl3]^ T . They carry quantum numbers B = Y = (1 mod 2)/3 and S = 1/2, and may 



even appear below the superconducting gap, a phenomenon unseen in normal superconduc- 
tors. In section 2, we will give arguments for the effective action in the CFL phase. In 
section 3, we show that superqualitons are topologically stable and identify their quantum 
numbers. In section 4, we give two complementary estimates for their mass and size. In 
section 5, we discuss their quantized spectrum. Our conclusions are in section 6. 



#1 We are suggesting that the notion of quark soliton or qualiton is more appropriate at high density than 
at zero density at which the qualiton picture has met with little success M 



2. In the CFL superconducting phase, quarks with opposite Fermi momenta pair 
with color and flavor locked. Model calculations pi| and effective theory calculations || 
have shown that this phase is favored energetically when the quark masses are light for 
a reasonable choice of parameters. The pertinent condensate will be taken in the (3, 3) 
color-flavor representation as 

ia Jb \ _ _/ ia jb \ _ ij abl 



1?A) = -{QRaQ J R ) = « e«e a »e a0I (1) 

where k is some constant, i,j are SL(2,C) indices, a, b are color indices, and a,f3 are 
flavor indices. By allowing for an arbitrary relative phase between the left- and right- 
condensate, parity could be spontaneously broken at high densities @. For finite k, both 
global color SU(3) C ^ and flavor SU(3)f^L,R symmetries are broken. The flavor-color 
locking in (|]) implies spontaneous breaking through the color-flavor diagonal. Specifi- 
cally: SU(3) C x (SU(3)l x SU(3)r) — > SU(3) c+ l+r, with the emergence of 8 pseudoscalar 
Nambu-Goldstone bosons together with 8 scalar Nambu-Goldstone bosons eaten up by glu- 
ons through the Higgs mechanism. There is an extra Nambu-Goldstone boson associated 
with U(1)b — ► Z% with no relevance to our discussion. 

To describe the low-energy dynamics of the color-flavor locking phase, we introduce a 
field Ul{x) which maps space-time to the coset space, Ml = SU(3) C x SU{3)l/SU{3) c+ l- 
Specifically, the left-handed field is 

1 7m 



ULaa(x) = lim ^eabcea^qf^-vp, x)q c I J j (v F , y), (2) 



where 7 m is the anomalous dimension of the diquark field of order a s and q(vp, x) denotes 
the quark field with momentum close to a Fermi momentum /jlvf [§] • The pairing involves 
quarks near the opposite edge of the Fermi surface. Similarly, we introduce a right-handed 
field Ur(x), also a map from space-time to Mr = SU(3) C x SU (3) r/ SU '(3) c +r, to describe 
the excitations of the right-handed diquark condensate. Under an SU(3) C x SU(3)l x SU(3)r 
transformation by unitary matrices (g c ,9L,9R), Ul transforms as Ul | -^ g*UL,g L and Ur 
transforms as Ul ^ 9*Ur9r- In the ground state of the CFL superconductor, Ul and Ur 
take the same constant value. QCD symmetries imply 

{ULaa} = ~ {URaa) = K5 aa . (3) 

The Nambu-Goldstone bosons are the low-lying excitations of the condensate, given as 
unitary matrices Ul{x) = g^(x)gL{x) and Ur(x) = g£{x)gii(x). For the present decom- 
position, we note the extra invariance under the (hidden) local transformation g c+ L+R{x) 
within the diagonal SU(3) c+ l+r through 

gj(x) -»■ gI(x)gl +L+R (x) gi,R( x ) -* gc+L+R(x)g LyR (x) (4) 



* 2 The breaking of global color is of course a misdemeanor that we will not address here. See [n5[ and 
Langfeld and Rho in M for discussions on this point. 



Hence, the spontaneous breaking of SU(3) C x (SU(3)l x SU(3)r) — ► 5C/(3) c+ l + r can be 
realized non-linearly through the use of Ul,r(x) or linearly through the use of g c ,L,R(x) 
with the addition of an octet vector gauge field transforming inhomogeneously under local 



9c+l+r{x)- This is the hidden local symmetry approach [21], in which the ensuing vector 
gauge field is composite and Higgsed, which is to be contrasted with a recent suggestion [22]. 
Clearly the composites carry color-flavor in the unbroken subgroup, with a mass of the order 
of the superconducting gap. Under some general conditions (e.g. vector dominance), the 



linear and nonlinear representations are the same [23|. Some of these points will be further 
discussed elsewhere. 

The current associated with det Ul (det Ur) is the left (right)-handed U(1)l (U(1)r) 
baryon number current. Because of the £7(1) axial anomaly, the field det Ul/ det Ur is 
massive due to instantons ^j. We decouple the massive field from the low energy effective 
action by imposing det Ul/ det Ur = 1. The Nambu-Goldstone boson associated with 
spontaneously broken U(1)b symmetry is described by det Ul ■ det Ur and responsible for 
baryon superfluidity. But since it is not directly relevant for our problem, we further choose 
det Ul ■ det Ur = 1 to isolate the Nambu-Goldstone bosons resulting from the spontaneous 
breaking of chiral symmetry, from the massive ones eaten up by the gluons. We now 
parameterize the unitary matrices Ul,r as 

U L (x) = exp (2in£ T A /F) , U R {x) = exp (2iUiT A /f) , (5) 

where T are SU{3) generators, normalized as Tr T T B = 8/2. The Nambu-Goldstone 
bosons II transform nonlinearly under SU(3) C x SU(3)l X SU(3)r but linearly under the 
unbroken symmetry group SU(3) c+ l+r- 

The effective Lagrangian for the CFL phase at asymptotic densities follows by in- 
tegrating out the 'hard' quark modes at the edge of the Fermi surface as suggested by 



Polchinski [24] and Weinberg [25]. This program has been carried out partly by one of 
us S. The effective Lagrangian for Ul,r is a standard non-linear sigma model in d=4 di- 
mensions and should include the interaction of Nambu-Goldstone bosons with colored but 



"screened" gluons G ^_. The SU(3) C current of Nambu-Goldstone bosons in the CFL phase 



consists of two pieces; one is from the Noether color-current and the other one is from the 
WZW term, given as 

JcL Afl = l ~F 2 Tr U L l T A d»U L + ^f^'Tr T A Ul x d v ULUl x d p ULUl x d a UL, (6) 

where the first term is the Noether current and the second one is from the WZW term. Ex- 
panding in powers of derivative, the effective Lagrangian for the (colored) Nambu-Goldstone 



* 3 The axial anomaly persists even in perturbation theory, and does not fully require the presence of 
instantons. The instantons of course persist even at asymptotic densities, albeit in a screened form. 

* 4 The magnetic gluons are Landau-damped in perturbation theory, but this may not be valid at asymptotic 
densities as suggested in [JL4|, if m* ~ (m^A) 1 ' 3 . 



bosons is then 

£= + ^fTi(8„ULd a ul)-^fli(8 i U L d,ul)+ a ,a-J cL +n L Cwzw 

In the presence of a chemical potential, only 0(3) is enforced in space. In carrying the 
derivative expansion, we assume the superconducting scale 4irA to be the large scale. (In 
effective theories the emergence of Air is purely a loop effect in d=4). The Wess-Zumino- 
Witten [pq , |27|| (WZW) term is nonlocal in d=4. Its normalization will be discussed below. 
The temporal and spatial decay constants Ft,s are fixed by the 'hard' modes at the Fermi 
surface. Their exact values are determined by the dynamics at the Fermi surface. For 
simplicity, we will assume F$ ~ Ft ~ F ~ A from now on. 

We note that (|^) only accounts for the gluon interaction to leading order in g s . To 
order g 2 s the gluon mediated interaction yields L, R mixing in a parity-invariant way. The 
breaking of chiral symmetry in the CFL phase via composite operators, e.g. {(qq) 2n ) ~ 
I {qq) | 2n > does not imply additional order parameters. No additional Goldstone modes are 
needed as all expectation values follow through simple Fierz rearrangements. 

The effective Lagrangian (0) in the CFL phase bears much in common with Skyrme- 
type Lagrangians [28, 29, |18| , with the screened color mediated interaction analogous to 



the exchange of massive vector mesons. Indeed in the CFL phase the 'screened' gluons and 
the 'Higgsed' gauge composites are the analogue of the massive vector mesons in the low 
density phase |l(| . In the CFL phase the WZW term is needed to enforce the correct flavor 
anomaly structure. Its structure is best described in 5-dimensional space as the homotopy of 
the Nambu- Goldstone manifold M, U^(M) = Z. Through the normalization Swzw = 2m7r 
with m £ Z, its coefficient ni and ur must be an integer to have a consistent quantum 
theory. 



3. Much like the effective Lagrangian for QCD at low density (giving rise to 
skyrmions), the low-energy effective Lagrangian in the CFL phase admits a stable (static) 
soliton solution, with a winding number given by the homotopy H^M) = Z . It is stable by 
the balance between the kinetic energy (attractive force) and the Coulomb energy (repulsive 
force) [pfl].^ Since the soliton quantum numbers are determined by the WZW term upon 



quantization, we need to determine the coefficient n. For skyrmions, Witten |26|] has shown 
that til and ur are fixed by the underlying flavor anomalies in QCD. Here we proceed to 
show that til and ur are fixed by analogous color-flavor anomalies in the CFL phase. 

First, consider the anomalous SU(3)l current, the pertinent U{\) subgroup of which 

* 5 There can be additional higher derivative terms that may significantly contribute to the stabilization 
of the qualiton which we shall not address here. Differently to the qualiton in the vacuum, one would also 
have to consider terms that are specific to the Fermi surface. 



is usually at the origin of the pion decay into two photons (i.e., ABJ triangle-anomaly). 
Since under the U(l) em electro-magnetic transformation the quark transforms as q i— > e ie ®q 
with Q = diag (2/3, —1/3. — 1/3), the £7(l) em transformation of the CFL condensate can be 
undone by a U(l)y hyper-charge transformation of SU(3) C with Y = —Q. Therefore, the 
unbroken U(1)q gauge boson (or modified photon) is a linear combination of the original 



photon and gluon [ 31 1 



Ap = A„ cos 6 + G Y sin0 (8) 

where cos 9 = g s /\/e 2 + g 2 and G Y is the gluon field for U(l)y C SU(3) C . Since the quark 
now couples to the modified photon with strength e = eg s / \Je 2 + g 2 , the unit charge of 
17(1)0 is e. 

A £7(l)/j transformation exp(ieQ) on U is nothing but a simultaneous transformation 
of g c = exp(iey) and g = exp (ieQ): 

U ^ e ieQ Ue- ieQ , (9) 

where we have used Y = —Q for the left multiplication. We see that the U(1)q transfor- 
mation Eq. (||) is identical to the C/(l) em transformation of the Nambu-Goldstone bosons 
associated with chiral symmetry breaking SU(3)lxSU(3)r i— ► SU(S)v in the QCD vacuum, 
except for the fact that the unit charge is now e = eg s /\/e 2 + g 2 instead of e. 

Therefore, all the modified electric charges of the Nambu-Goldstone bosons in the 
CFL phase are integral multiple of e, the modified electric charge of the electron. Note 
that all the quarks and gluons have integer (modified) electric charges too, since quarks and 
gluons transform as Eq. (||) under U{1)q |l5|]. Specifically, 

Gf, ^ e^G^e-^ = e ieY G IM e~ leY (10) 

q ^ e i€ Qqe ie ® = e itY qe itQ (11) 

The SU(3)l anomaly is now given as 

~2 

A_ e r TX ,fr r A/S2 



d^4 = ^2 TV {T A Q 2 ) e^F^FO", (12) 

where F^ v is the field strength tensor of the modified photon associated with U(1)q. Note 
that there is no color or flavor factor in the coefficient because neither color nor flavor 
is running around the triangle diagram responsible for the SU(3)l anomaly. In the CFL 
phase color and flavor are locked, and the modified pion and photon carry color. In the 
CFL phase, 

/0 -1 -l\ 



Qc 



(13) 



1 

Vi o o j 

and only three quarks contribute to the anomaly in the Jjf current: Tr (T 3 Q 2 ) = 1/2 
(l 2 + l 2 ) + (-1/2) (-1) 2 = 1/2. 



To fix the coefficient n of the WZW term from the SU(3)l anomaly, we introduce 
the modified electro-magnetic interaction by gauging the U(1)q symmetry in the effective 
action (|?]) through 

d»U » D»U aa = d^Uaa - ieA^ (Y ab U ba + U af3 Qp a ) . (14) 

As shown by Witten pq], the U{\)q invariant WZW term then contains a term 

S W zw(U, A,) B -^rfe^F^F^, (15) 

which agrees with (12) if m, = 1. Similarly, one can also show that ur = 1. 

On the other hand, the topologically conserved current for the soliton in (|?]) is 

Jl = ^2^ upaTT ui l d u u L uz l d p u L ui l d a u L . (ie) 

Since the sigma model description of the SU(3)l quark current J A = qLT A ^^qL contains 
an anomalous piece from the WZW term 

J t ^ ^^"^ T A V I r x d v U L U I r X d p V L V L - x d <T U L (17) 



and the topological current Eq. fllq) corresponds to the anomalous piece of the U{\)l current 
in the a model description p6| , |32] , |33], we find the soliton of unit winding number has a 
quark number 1 or baryon number B = 1/3. This is Kaplan's qualiton |l^] in the CFL 
phase. We are referring to it as the superqualiton to distinguish it from Kaplan's qualiton in 
the vacuum. Since in the CFL phase baryon number is spontaneously broken £/#(l) — > ^2> 
it is clear that B = 1/3 is the same as B = 1, hence B = (1 mod 2)/3. The fact that til = 1 
(or riR = 1) implies that the superqualiton is a fermion. The superqualiton carries positive 
parity, since the CFL condensate is parity even. Further quantum numbers along with the 
excited spectrum of the superqualiton will be discussed below. 



4. Some important questions regarding superqualitons concern their mass Ms and 
their interactions. Indeed, if M$ happens to be less than A (where A is the superconductiv- 
ity gap), then superqualitons may exist inside the superconducting CFL gap with ordinary 
baryon quantum numbers, a phenomenon unseen in conventional superconductors. Attrac- 
tive interactions among light superqualitons may even trigger a total rearrangement inside 
the QCD superconductor. 

A simple estimate of the mass follows from scale (virial-like) arguments that are in 
general not specific to the detailed form of (j^). Indeed, for stable static configurations a 
typical gradient in (|7|) is d ~ A. Since F ~ A and the size of the superqualiton Rs ~ 
27r/(4vrA), then 

M 5 ~Qf 2 A 2 ) (!*!§) ~^A (18) 



which is less than A/2. Both the mass and the size are fixed by A. To create 2 superquali- 
tons costs less than to excite 2 bare quarks from the CFL phase. Such superqualitons 
may either condense (crystallize) into a new superconductor with a reduced gap, or exhibit 
a liquid or gas phase with either normal or superconducting properties. In order to ad- 
dress these issues, the specific dynamics of the many-superqualiton interaction is required. 
Unfortunately we have presently no good understanding of this aspect of the problem. 

A more conservative estimate of the mass of the superqualiton can be made by relying 
on the effective Lagrangian (0). Following Skyrme |28| and others |ll| (and references 
therein) we seek a static configuration for the field Ul in SU(3) by embedding an SU(2) 
hedgehog in color-flavor in SU(3), with 



U Lc {x) 



o if-r0(r) 



, U R = 0, Gl = u(r), all other G~ = 0, 



(19) 



where r's are Pauli matrices. The radial function 9(r) is monotonous and satisfies 

6(0) = 7T, 0(oo) = (20) 

for a soliton of winding number one. (Note that we can also look for a right-handed soliton 
by switching off the Ul field. The solution should be identical because (^) is invariant under 
parity.) This configuration has only non- vanishing color charge in the Y direction 



J, 



Y 



sin 



2vrr 2 



(21) 



and all others are zero. As shown in [BO] , the energy of the configuration (fi^) is given as 



E[ui,i 



Airr dr 



1/2 9 / ,i smf 

--J 2 + F 2 [e ' 2 + 2— n 



The total charge within a radius r is 



\Y , 



tY,J\ 



j^ij 



Q Y ( r ) =g 8 Tr YJJ (r f )4nrr f dr' 
Jo 



-9s 



9s V . 2 

+ 2^ Sm ' 



(r) — sin 9(r) cos 6(r) — it 



TV 



Using Gauss law with a screened charge density, we can trade uj in terms of 9(r) 

Q Y (r) 



UJ 



-rtiET 



Airr 2 



(22) 



(23) 



(24) 



where rriE = y/6a s /ir fj, is the electric screening mass for the gluons (note that the magnetic 
gluons are not needed for the static configuration). Hence, the energy functional simplifies 
to 



E[9] 



£(r) dr = / 47rdr 
o Jo 



F 2 ( 



r 2 9' 2 + 2 sin 2 . 



+ 



Ots_ 

2vr 



8111 V COS ( 



7T 



2r 



-2mET 



(25) 



where a s = <7 2 /(47r). The squared size of the superqualiton is R s = (r 2 ) where the averaging 
is made using the (weight) density £(r). The equation of motion for the superqualiton profile 
9{r) is 

(r 2 9'j = sin29 + -— sin 2 g tt - 9 + -sin2fl (26) 

subject to the boundary conditions (pGJ). 

We have solved ( |26|) numerically for several values of m_E and a s . Given a fixed 
a s = 1, we find the soliton mass Ms = 1.53, 1.06, 0.75 x 4ttF and the soliton radius 
Rs = 0.29, 0.20, 0.15 F _1 for m^/F = 2,5,10, respectively. For various values of a s = 
1.0, 0.5, 0.2 with a given m E /F = 2, we get M s = 1.53, 1.21, 0.87 x 4vrF and R s = 
0.29, 0.24, 0.17 F . We see that both the mass and the size of the soliton decrease as 
the repulsive Coulomb force decreases, showing that the soliton gets smaller for the weaker 
coupling, since the color-electric force, which balances the kinetic energy of the soliton, is 



less repulsive [30|. For the high density region, where the QCD coupling is small, we expect 
the higher order derivative terms to become more important in stabilizing the soliton, since 
the Coulomb forces become negligible following screening (an alternative would be through 
the use of gauge composites). This trend is expected from the screened character of the 
interaction, and seems to differ from the one we suggested earlier by our simple scaling 
arguments. Which scenario (among the various alternatives) is favored will depend on the 
content of the effective Lagrangian implied by first principles of QCD which may not be 
reliably embodied in (0). 



5. To access the quantum numbers and the spectrum of the superqualiton, we note 
as usual that for any static solution to the equations of motion, one can generate another 
solution by a rigid SU(S) rotation, 

U(x) i-» AU(x)A~ 1 , A G SU{3). (27) 

The matrix A corresponds to the zero modes of the superqualiton. Note that two SU(3) 
matrices are equivalent if they differ by a matrix h € U(l) C SU(3) that commutes with 
SU(2) generated by t ® I. The Lagrangian for the zero modes is given by substituting 
U(x,t) = AtyUcfflAfy)- 1 61. Hence, 



L[A] = -M s + -/ Q/3 Tr r at A- 1 iTt T p A' x A - i-Tt YA~ X A, (28) 

where I a p is an invariant tensor on M. = SU(3)/U(1) and the hyper charge Y is 

/l \ 
Y = ~ 1 . (29) 

\0 -2/ 



Under the transformation A(t) i— > A(t)h(t) with h £ U(1)y 

L i-» L - -Tr YhT x h. (30) 

Therefore, if we rotate adiabatically the soliton by 9 in the hypercharge space in SU(3), 
h = exp(iY6), for time T — > oo, then the wave function of the soliton changes by a phase 
in the semiclassical limit: 

7p(T) ~ e^^VCO) = e ie/3 V>(0), (31) 

where we neglected the irrelevant phase —M$T due to the rest mass energy. The simplest 
and lowest energy configuration that satisfies Eq. ( |3l| ) is the fundamental representation of 
517(3). Similarly, under a spatial (adiabatic) rotation by around the z axis, h = exp(ir 3 #), 
the phase of the wave function changes by 9/2. Therefore, the ground state of the soliton 
is a spin-half particle transforming under the fundamental representation of both the flavor 
and the color group, which leads us to conclude that the soliton is a massive left-handed 
(or right-handed) quark in the CFL phase. 

It is interesting to note that the phase continuity between high-density quark matter 
and low-density nuclear matter, conjectured by Schafer and Wilczek |16|], carries to the 
Skyrme picture of baryons. One way to address the specifics of the continuity besides a 
spectrum analysis, is via anomaly-matching. This issue will be addressed elsewhere. 



6. In the CFL phase of the QCD superconductors, baryons emerge as qualitons 
(superqualitons). We have shown that the quantum numbers of the superqualitons are 
dictated by a WZW term suitably normalized by the CFL triangle-anomaly We have 
presented general arguments for the effective Lagrangian in the CFL phase. Superqualitons 
can be stabilized by screened colored electric interactions for a variety of chemical potentials, 
exhibiting a rich and colorful spectrum. Our estimate of the mass is on the higher side 
of the superconducting gap, although simple scaling arguments show it otherwise. Our 
calculations do not exhaust all parameter space for effective model calculations in QCD 
superconductors, and point at the necessity of a first principle assessment of the effective 
Lagrangian (0). One fascinating possibility suggested by the superqualiton description of the 
hadron-quark continuity is that light fermionic modes with M s < A get excited within the 
CFL superconducting gap, a phenomenon unseen in normal superconductors. On general 
grounds, the interaction between light superqualitons may cause a rearrangement in the 
superconducting phase with important implications for the bulk properties such as the 
specific heat. The latter is paramount to cooling curves of neutron stars. 
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Note added 

After posting our paper, |jj| appeared in which some of the present issues were 
discussed. We note that the linear interpretation of the symmetry breaking pattern adopted 



by the authors in [35] differs from ours and that our nonlinear description does contain L, R 



mixing to order gj contrary to their statement. 
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